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Abstract 

A two-electron system confined in two coupled semiconductor quantum dots is investigated as a 
candidate for performing quantum logic operations on spin qubits. We study different processes of 
swapping the electron spins by controlled switching on/off the exchange interaction. The resulting 
spin swap corresponds to an elementary operation in quantum information processing. We perform 
a direct time evolution simulations of the time-dependent Schrodinger equation. Our results show 
that - in order to obtain the full interchange of spins - the exchange interaction should change 
smoothly in time. The presence of jumps and spikes in the corresponding time characteristics leads 
to a considerable increase of the spin swap time. We propose several mechanisms to modify the 
exchange interaction by changing the confinement potential profile and discuss their advantages 
and disadvantages. 
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I. INTRODUCTION 



In quantum computation, one and two-qubit logic gates play a crucial role, since they are 
allow us to perform arbitrary quantum logic operation Recently, a practical realization 
of these gates is a challenge for many physical laboratories. Among several propositions 
of constructing the systems performing the two-qubit gates, the most promising are these 
based on the mechanism of controlled switching on /off the exchange interaction between 
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qubits q, y, y, y, y, y, y, y, q, 
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14j . Tuning the exchange interaction 



between electrons in coupled quantum subsystems can lead to the interchange of qubits and 
the performance of a designed quantum logic operation. 

A physical implementation of the controlled exchange interaction is itself a difficult task. 
Moreover, the elements of a future quantum computer have to fulfill several conditions, 
known in the form of DiVincenzo criteria 15] . One of them is a scalability, which allows 



to join the elements into a large computational machine. It is expected that semiconductor 
nanodevices, in particular those consisting semiconductor quantum dots (QDs), should be 
scalable to large enough size. Recently, the electron spin states are the most promising 
candidates for qubits 6|, |8|, |9j due to the long coherence time 16j . The quantum information 
processing can be performed via the changes of electron spin states js), The spin of 
the electron can be rotated as a result of precession in the external static magnetic field 
0, [3] . The spin swapping can also be obtained when irradiating the electron system by 
the electromagnetic wave with the frequency adjusted to the Zeeman splitting in the external 
magnetic field jsl. 

In the present paper, we study the process of swapping the electron spins, which results 
from the controlled switching of the exchange interaction between the electrons in coupled 



ms m coupled 

)inQ,y,laQ, 



QDs. The tuned exchange interaction should lead to the rotation of electron spin 
[l7,0|. The theoretical model proposed takes into account all the three space dimensions and 
the electron spins. We apply the adiabatic approximation in order to decouple the transverse 
and longitudinal degrees of freedom and perform the integration over the transverse spatial 
coordinates. This leads to the effectively one-dimensional two-electron problem [2l|, which 
can be solved by the numerical method with an arbitrary precision 22]. The system under 
study is described by the two-particle wave function of the form of four-component state 
vector, which takes into account all possible spin configurations. We consider different 



2 



methods of turning on and off the exchange interaction between the electrons and investigate 
the resulting rotation of electron spins. The present computer simulations are based the 
accurate numerical solutions of the time-dependent Schrodinger equation, i.e., they allow us 
to trace the time evolution of the electron system in a direct manner. Therefore, the present 
results fill in a gap between the quantum information theory and experiment. We hope that 
these results will serve as a guide for the experimental groups who are involved in designing 
and constructing the nanodevices with spin qubits. 

The paper is organized as follows: the theoretical model is presented in Sec. II, the 
numerical method is described in Sec. Ill, the results are presented in Sec. IV (for the 
vertically coupled QDs) and in Sec. V (for the laterally coupled QDs). Section VI contains 
conclusions and summary. 



II. THEORETICAL MODEL 



The model applied in the present pa per is an extension of that proposed by the present 



authors in Ref. 



23| . Contrary to Ref. [23j], we explicitly take into account the spin states 



of electrons. For sake of completeness, we repeat below the major steps of the approach 



231 ] . We study the two electrons localized in the double QD nanostructure. We assume 



the cylindrical symmetry of the system with symmetry axis (z axis) going through the 
geometrical centers of the QDs. The potential confining the electrons in the x — y plane 
is taken to be so sufficiently strong that the differences between the energy levels resulting 
from the x — y space quantization are much larger than the energy of spin swapping. This 
potential is usually called the lateral confinement potential. We approximate the lateral 
potential by the two-dimensional harmonic oscillator potential and assume that both the 
electrons occupy the ground state in the x — y motion. In other words, we assume that the 
lateral electron degrees of freedom are frozen. The above assumptions allow us to reduce the 



starting three-dimensional problem to the effectively one-dimensional problem [21( with the 
electron-electron interaction being the following function of z x and Z2 coordinates of both 
the electrons: 

U eff ( Zl - z 2 ) = ^^e^--l 2 erfc (yfp\z x - z 2 \) , (1) 

where e is the vacuum electric permittivity, e is the static electric permittivity, (3 = 
m e hu)±/2, me is the electron effective conduction band mass, and hu)± is the excitation 
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energ y le vels of the electron transverse motion. The Hamiltonian of the system takes the 
form 1 2 11 



H 



h 2 



2m P 



^2 



(2) 



C)7 2 

where U(zi) is the potential energy of the vertical confinement of the ith electron. In the 
following, vertical confinement U(z) will be taken in a form of two potential wells separated 
by the potential barrier. The energy is measured with respect to the conduction band 
bottom of the quantum well material. In the calculations, we take on the parameters, which 
correspond to nanostructure based on GaAs: m e = 0.067m e0 , where m e0 is the free electron 
rest mass, and the static electric permittivity e = 11.0. 

In the two-electron system, we are dealing with four independent spin states; therefore, 
the total wave function of the system can be represented by the following four- component 
vector: 

where <x,- (j = 1,2) are the spin coordinates of the electrons, ip nm = ^ n m{z\-,z^) are the 
basis wave functions with the indices n and m — f and j, which correspond to the z-spin 
component eigenvalues +H/2 and —H/2. Basis wave functions ip nm do not possess a well 
defined symmetry. However, total wave function ([3]) has to be antisymmetric with respect 
to simultaneous exchange of the space and spin coordinates. This property imposes the 
following conditions on the basis wave functions: 



(3) 



i/>ll(zi,z 2 ) = -ip n (z 2 ,zi) , 
^li(zi,z 2 ) = -^ u (z 2 ,zx) . 



(4a) 
(4b) 
(4c) 



In representation ([3]), the operators of the total spin are expressed by the following 4x4 
matrices: 
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where 1 denotes the unit 2x2 matrix and a x ,a y ,a z are the spin Pauli matrices. The cor- 
responding expectation values are calculated as (cr£°*) = (^o^ *^), where k = x,y,z. Mul- 
tiplying (cr£) by h/2 we obtain the expectation values of the components of the total spin, 
i.e., S£* = (h/2)(oj«). 

Due to their indistinguishability the electrons can not be numerated; therefore, we can 
not determine which electron is in a given spin state. However, we can distinguish the two 
quantum dots, determined by the corresponding potential wells. We shall call them the 
left (L) and right (R) quantum dot (potential well). If both the QDs are separated by the 
potential barrier, which prohibits a tunneling of electrons, then in the ground state of the 
system, each electron is localized in a single QD. Owing to this, we can determine the spins of 
electrons in the left and right QD. For this purpose, we introduce the auxiliary wave function 
(p(zi, Z2), which will be called the reference wave function. This wave function is a solution 
of the Schroinger equation, for the two distinguishable particles, which - with an exception 
of indistinguishability - possess all the properties of the electrons. For the distinguishable 
particles we do no perform a symmetrization of the two-particle wave function. Instead 
the solutions found for ip[zi,Z2) correspond to the configurations, in which one electron 
(described by coordinate z{) is localized in the left QD and the other (z 2 ) is in the right 
QD. This state always exists if the tunneling through the barrier is not possible. In each 
simulation performed, these configurations are realized at the initial (t = 0) and final (t = T) 
time instants. The reference wave function will serve to a construction of the initial state 
wave function and a determination of the spin states of electrons in both the QDs. The 
expectation values of spin components are calculated as follows 



h h 
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where j = L,R corresponds to the electron localized in the left (L) and right (R) QD, and 

\(<p\MJ 

The spin matrix operators in Eq. (jBJ) have the form 

l\ r [0 -l\ T ( 1 
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The operator of the fcth component of the total spin is defined as Sl ot = (h/2)(a^ + aj}). 

III. SETTING UP THE COMPUTER EXPERIMENT 

We simulate the process of swapping the spins for two electrons localized in two QDs, 
which are separated by the potential barrier. We consider both the vertical and lateral 
geometry of the QD nanostructure. In the process studied, the spin of each electron localized 
in the left or right QD is flipped. The spins are swapped as a result of the controlled switching 
of the exchange interaction between the electrons. We discuss two methods to control the 
exchange interaction. According to the first method, the exchange interaction is switched 
on (off) by lowering (rising) the potential barrier. The second method is based on the 
changes of the depth of one of the potential wells, which leads to the flow of electrons and 
their localization in the same QD. In the QD-based nanodevices, both the methods can be 
implemented by changing the external voltages, which are the source of the electrostatic 
field forming the coupled QDs [3], or by locating the QD system in the electromagnetic 
field J 

Based on the analogy with Rabi oscillations in a two-level quantum system, the process 
of switching on/off the exchange interaction that leads to the change of the initial spin 
orientation to the opposite one, will be called the 7r pulse. This process lasts for time 
T n . We shall also deal with processes, after which the expectation value of z electron spin 
component does not reach ±h/2, i.e., < I n such processes, the change of spin 

orientation is not complete. In electrostatic QDs, the 7r pulse can be realized by the proper 
change of gate voltages. 

The additional purpose of the present study is the optimization of nanodevice parameters 
and time changes of the confinement potential in order to make the spin interchange time 
possibly short. The minimization of T w is important since this time interval determines the 
duration of elementary quantum logic gate operation that should be much shorter than the 



coherence time of qubits. It is required [15] that the ratio of the gate operation time to the 



coherence time should be less than 10 4 in order to complete a computation and several 



* i = V(z 1 ,<Ti,Z2,<T 2 ,t = 0) 



(9) 



error correction runs [24| before the decoherence destroys the information stored in qubits. 
In GaAs, the spin coherence time has been estimated to be longer than 1 ms j^J. 

At the initial time instant in each simulated process, the system is always in the singlet 
ground state, in which the electron localized in the left (right) QD possesses the spin z 
component S%(t = 0) = +h/2 (S*(t = 0) = -h/2). The two-electron wave function that 
satisfies these initial conditions has a form 

/ \ 

V o / 

where (p(zi,z 2 ) is the reference wave function defined in Sec. II. Wave function (Q fulfills 
symmetry constraints (j3J). After a normalization, the wave function (jUJ) is subjected to a 
temporal evolution according to the time-dependent Schodinger equation 

d^(z u z 2 ,t) 

in ^ = HV(zi,Z2,t) , (10) 

where H is given by Eq. (j2j) and the time dependence of potentials U(zi,z 2 ,t) and 
U e ff(zi, z 2 , t) will be determined in the following sections. 

The spin expectation values are calculated using Eq. ([6]) as follows. First, the reference 
wave function (p(zi, z 2 ) is calculated by the imaginary time step method [22|. The reference 
wave function is next applied to construct total wave function (Q for t = 0. Wave 
function is used as the initial condition to perform the time simulation of Schrodinger 
equation ( TlOl) . In the last step, we apply (p(zi,z 2 ) to calculate u (Eq. ([Zj) and S k (Eq. (jHJ). 

The results of simulations are presented as a function of duration time T of the process 
of switching the exchange interaction. We note that we are dealing with two time intervals 
T and T 71 ", which have a different meaning. Time interval T determines the duration of 
the arbitrary process of changing the exchange interaction. Time interval T n determines 
the duration of such process, during which the spins are fully swapped. At the end of 
each simulation, we record the expectation values of z spin component of the electron in 
the left (L) and right (R) QD. These final values of spins are denoted by Sf' R , where 
Sj' R = S L ' R (t = T). Moreover, we record the energy difference AEfi = Ef — Ei between 
the energies of the final {Ef) and initial [Ei) states of the system. The non-zero value of 
AEfi = Ef — Ei gives the evidence that the process in non-adiabatic. 



The numerical procedure applied allows us to observe the spins of single electrons in an- 
other manner. Namely, we can directly determine the temporal evolution of the second (/0TJ.) 
and third (ipii) component of the total wave function ([3]), whereas - due to the symmetry - 
it is sufficient to detect one of these components. The full swapping of spins corresponds to 
the interchange of these two components, i.e., 



V o / 



V o / 



(ii) 



IV. VERTICALLY COUPLED QUANTUM DOTS 



The vertically coupled QDs |26| can be described by the simple model with two rectan- 
gular potential wells separated by the rectangular barrier (Fig. 1). The vertical QDs were 



271 ] . The exchange 



fabricated on the basis of GaAs/InGaAs/AlGaAs heterostructure 
interaction is switched on by lowering the height of potential barrier to the potential well 
bottom. The exchange interaction is switched off by raising the potential barrier to the 
initial height. The time dependence of the process of changing the barrier height can chosen 
in many ways. In the present paper, we have chosen the time dependencies described by the 
step, linear, and smooth functions. 

We choose such values of the parameters of the nanostructure (Fig. 1) that the electron 
tunneling is negligibly small in the initial and final time instant. Then the electrons are 
well separated in space and localized in the left (right) QD. Owing to this the spin states 
of both electrons can be exactly known in both the initial and final instant. We choose the 
same widths of the potential wells and potential barrier cIl — d,R — ds — 10 nm and the 
potential well depth Vl = Vr = —150 meV. At the initial time instant, we take on the top 
barrier energy Vb = (the conduction band bottom of the barrier material is chosen as the 
reference energy and set equal to zero). 
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FIG. 1: Confinement potential profile E = V as a function of vertical coordinate z and energy Vb 
of the top of potential barrier, (a) Vb = 0, (b) Vb = —50 meV, (c) Vb = —100 meV, and (d) 
V B = -150 meV. 




A. Time changes of the potential barrier 

First, we study the effect of the time step dependence of the barrier height [Fig. 2(a)]. In 
this process, the top of the barrier is rapidly lowered from Vb — to Vb — —150 meV, i.e., the 
potential well bottom. We let the system to evolve for time t = T, after which the barrier is 
rapidly raised to the starting value. We have performed a series of simulations for different 
operation time T. At the beginning and the end of each process, we have determined z spin 
component expectation values S^' f R . 



(a) (b) (c) 




t [ps] t [ps] t [ps] 

FIG. 2: Time dependence of top barrier E = Vb changes: (a) step function, (b) linear function, 
(c) smooth (cosinus) function. 

The results, displayed in in Fig. 3(a), show that the spin expectation value detected at 
the final instant of each process rapidly oscillates as a function of time duration T of the 
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process. However, the amplitude of these oscillations does not reach k/2. Therefore, the 
sudden changes of the barrier do not lead to to the full interchange of spins at any time 
interval T studied. The behaviour of energy difference AEfi between the final and initial 
state [Fig. 3(b)] allows us to explain this effect. Similar to the average electron spin, energy 
difference AEfi exhibits rapid oscillations as a function of duration time T of the process. 
We interpret the results of Fig. 3(b) as follows: during the sudden jumps of the barrier 
height, the electron system makes transitions to the excited states. As a consequence, the 
energy of electrons in the final state increases, i.e., the energy difference is always larger than 
zero and exhibits jumps, and the electron wave function spreads out over the two potential 
wells. Therefore, the switching of the interaction between the electrons in a rapid manner 
is not an adiabatic process. 



(a) (b) 




T [ps] T [ps] 

LP L R 

FIG. 3: (a) Expectation value Sy = Sf' of z spin component detected at the end of the process 
of switching on and off the exchange interaction as a function of its duration time T for step-like 
changes of the barrier. Solid (dashed) curve shows the results for the electron in the right (left) QD. 
(b) Energy difference AEfi between the final and initial states energies as a function of duration 
time T of the process. 



Next, we study the two-step process, during which the potential barrier is a linear function 
of time [Fig. 2(b)]. In the first step, for < t < T/2, the top barrier energy is lowered 
according to 

V- /jA \rmax c\(\rmax \rmin\ ^ /10\ 
B{t) = V B - 2,{V B - V B )— , (12) 

where V B ax (V B tn ) denotes the largest (smallest) energy of the barrier top in time interval 
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T. In the second step, for T/2 < t < T, the potential barrier is raised according to 

V B (t) = V™ n + 2(Vr x - V™ n ) - . (13) 

The results (Fig. 4) show that spin expectation values S^' R are periodic functions of time T. 
Each maximal (minimal) value S R = +H/2 (Sj = —h/2) of z spin component corresponds 
to the state, in which both the electrons possess the spins of the opposite orientation with 
respect to that in the initial state. Therefore, there exists a series of duration times T, 
for which the z spin component of each electron changes its sign, i.e., the electrons fully 
exchange their spins. The subsequent extrema in Fig. 4(a) correspond to the 7r pulses with 
duration time T£ = nT^ + const, where index n numerates the extrema. If we increase 
process duration time T, we simultaneously increase the time interval during which both 
the electrons are localized in the same region of space. As a consequence, for the sufficiently 
long time T the spin swapping can occur many times. If the process lasted for infinitely 
long time, the spin interchange would occur infinitely many times. This scenario would 
be realized for an isolated system (non-interacting with an environment). In this case, 
spin expectation values Sf' R would oscillate with amplitude h/2 infinitely many times, i.e., 
there exists an infinite series of time intervals T£, for which the spins are fully interchanged. 
However, in real systems, we deal with the decay of quantum states, which leads to the energy 
dissipation, and the decoherence, which randomly changes the relative phase of spin qubits. 
These processes will lead to the decreasing amplitude of S^' R for increasing time interval 
T. Therefore, in order to perform successful quantum logic operation the spin swapping 
process should last for possibly short time. However, there exists the lower limit on time 
interval T, during which the barrier is changed. This limit results from the requirement of 
the adiabaticity of the process: the barrier has to be changed sufficiently slowly, i.e., in an 
adiabatic manner, in order to leave the system in the ground state. If time interval T is too 
short, the system can be excited to higher-energy excited states, which leads to the increase 
of its energy and the non-adiabaticity of the process. 

The non-adiabaticity of the spin interchange process can be also observed in the behaviour 
of Sj ,R (T) [Fig. 4(a)]. A closer look at the present computational results has led us to the 
observation that the first maximum is slightly lower than the other maxima. For the first 
maximum the exact spin expectation value reaches h/2 with accuracy 98.7%. Similarly, due 
to the symmetry, for the first minimum, Sj reaches the value — h/2 with the same accuracy. 
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The extrema, which correspond to the second % pulse, reach ±h/2 with accuracy 99.4%. The 
next extrema tend to ±h/2 with a high accuracy. Therefore, we deal with several processes 
of incomplete interchange of spins before Sy start to switch between ±h/2. The lowering 
of the amplitude of Sj' R is more pronounced if the process of changing the potential barrier 
is more rapid. The present results show that several processes of incomplete spin swapping 
can occur before we achieve the full interchange of spins, which leads to an elongation 
of the spin exchange time. We note that this effect is not expected within the model, 
based on the effective Heisenberg Hamiltonian. In this model, the exchange interaction 
between the electrons is described by the Hamiltonian H(t) = J(t)Si ■ S 2 , where J(t) is 
exchange interaction energy and the spin operators Sj (j = 1,2) are expressed in terms 
of the corresponding Pauli matrices. Our estimates show that - in real quantum systems 
- the spin interchange time will be longer than that resulting from the Heisenberg model, 
which do not take into account the physical implementation of switching on/off the exchange 
interaction. 



(a) (b) 
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FIG. 4: (a) Expectation value ,P = S^ ,R of z spin component detected at the end of the process 
of switching on and off the exchange interaction as a function of duration time T for linear changes 
of the barrier. Solid (dashed) curve shows the results for the electron in the right (left) QD. (b) 
Energy difference AEfi between the final and initial state energies as a function of duration time 
T of the process. Inset zooms in the part of this figure. 

Fig. 4(b) displays the dependence of energy separation AEfi on time interval T. In the 
inset, we zoom in this dependence for the narrower scale of T, which allows us to observe the 
oscillations of AEfi. The amplitude of these oscillations decreases if the duration time of 
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spin interchange increases. Moreover, AEfi asymptotically approaches zero for large T. For 
small T the positions of the local minima of AEfi do not coincide with T£ [cf. Fig. 4(a)]. 
For the chosen nanostructure parameters, this coincidence appears for the third tt pulse at 
T = 8.2 ps. For linear changes of the barrier we deal with rapid turning on the process at 
T/2, which causes the non-adiabaticity of this method of spin interchange. 

Based on the above results, we expect that we will obtain better results if the temporal 
changes of the potential barrier will be modelled by a smooth function of time. In the 
following, we study the smooth time dependence of the potential barrier, which is modelled 
by the following function [Fig. 2(c)]: 

V B (t) = l(Vr x ~ V™ n )(l - cosut) + , (14) 

where u = 2n/T. The results of simulations [Fig. 5(a)] show that the amplitude of S^ ,R (T) 
approaches h/2 much faster than for step-like and linear changes of the barrier. Already 
the second maximum reaches h/2 with accuracy 99.99%. Moreover, based on the results for 
energy difference AEfi [Fig- 5(b)] we can regard this process to be adiabatic for T > 2 ps. 
The very fine jumps of AEfi near zero, shown in the inset of Fig. 5(b), result from the small 
numerical errors. Therefore, the amplitude of the periodic changes of Sf ' R remains constant 
(and equal to h/2) for T > 2 ps. 

The spin of the electron localized in the left (right) QD can be determined from Eq. (J6j) 
only if the probability of finding both the electrons in the same QD is zero. This state is 
realized at the initial and final time instant in each process of the barrier changes. At the 
intermediate time instants, the electrons are not spatially separated; therefore, quantities 
S L,R = Sz' R , calculated from Eq. (jHJ), are not equal to the eigenvalues of z spin component 
of the electron in the left and right QD. Nevertheless, for arbitrary time, S L ' R can serve as 
a control spin index. 

In Fig. 6 we present the time dependence of control spin indices and energy of the system, 
which correspond to the pulses with duration times T[ and T£ for linear and smooth changes 
of the barrier. The energy is measured relative the ground state energy Eq of two electrons 
in the lateral confinement potential (Eq = 2hu±, where hu± = 40 meV). The behaviour of 
these time characteristics is different for the smooth and non-smooth changes of the barrier. 
In the case of the non-smooth process, the discontinuity of the first time derivative of the top 
barrier energy for t = T/2 leads to the corresponding discontinuity of the energy versus time 
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FIG. 5: (a) Expectation value S^' F = Sj ' R of z spin component detected at the end of the process 
of switching on and off the exchange interaction as a function of duration time T for smooth changes 
of the barrier. Solid (dashed) curve shows the results for the electron in the right (left) QD. (b) 
Energy difference AE/i between the final and initial state energies as a function of duration time 
T of the process. Inset zooms in the part of this figure. 



plot [cf. dotted curves in Figs. 6(a) and 6(b)]. The consequences of this discontinuity can 
also be seen in the time dependence of the control spin indices [cf. solid and dashed curves 
in Figs. 6(a) and 6(b)]. For the pulse with time duration T[ = 1.58 ps the spin indices 
begin to oscillate after reversing the time evolution of the potential barrier [Fig. 6(a)]. The 
amplitude of these oscillations decreases with the increasing barrier height and falls down 
to zero when the barrier height reaches the starting level. For the three times longer pulse 
T£ = 4.91 ps reversing the time changes of the barrier leads to smaller oscillations [cf. Fig. 
6(b)]. Nevertheless, the time dependence of the control spin indices is still slightly perturbed. 

In the case of smooth time changes of the barrier [Eq. (|14p ]. the energy of the system 
is also a smooth function of time [Figs. 6(c,d)]. In spite of this, for the first 7r pulse with 
Ty = 0.99 ps the changes of the potential barrier are still seen by the system as rapid changes, 
which leads to bending the spin characteristics [Fig. 6(c)]. However, the time duration of 
the second n pulse (T^ = 3.04 ps) is sufficiently long so that the spin characteristics [Fig. 
6(d)] stay smooth during the full process of switching on/off the exchange interaction. In the 
cases of linear and smooth changes of the potential barrier, the time evolution of spin indices 
S L ' R (t) for the second 7r pulse with time duration T£ suggests that the full interchange of 
spins occurs as a result of many-fold exchange of spins between the electrons. For the 
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FIG. 6: Spin control indices S^' P = S L,H and energy E of the two-electron system as functions of 
time t. Plots (a) and (b) [(c) and (d)] correspond to the linear [cosinusoidal] changes of the potential 
barrier, plots (a) and (c) [(b) and (d)] correspond to the first [second] ir pulse with duration time 
T[ [Trf]- Solid (dashed) curves display the control spin index of the electron in the right (left) QD 
and the dotted curves show the energy of the system. 



second pulse the curves S L ' R (t) go through zero three times [cf. Figs. 6(b,d)]. Thus, we 
conclude that in these processes we are dealing with the three-fold incomplete swapping of 
spins before we finally obtain the full interchange of spin orientations. In order to obtain 
an additional support for the occurrence of this process, we have observed the time changes 
of the second and third components of the total wave function [Eq. (j3j)]. The results of 
the simulations performed show the many-fold interchange of the second and third wave 
function components, which means that - in fact - the spins of the electrons change many 
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times. 

Based on these results, we have chosen to further studies the method of the smooth 
switching of the exchange interaction, which seems to be the most promissing in a physical 
realization. 



B. Effect of quantum dot size and potential well asymmetry 

If the size of the coupled QD nanostructure is sufficiently large, the Coulomb interelectron 
repulsion leads to a localization of electrons at the opposite potential well boundaries [32] . 
Then, the electrons form the Wigner molecule 32| . In the case of Wigner localization, the 



lowering of the potential barrier can cause no interchange of spins during a time, which is 
shorter than the spin coherence time. We have checked the possibility of the occurrence of 
this effect by investigating the influence of the QD size of the time duration of spin inter- 
change process. The calculations have been performed for the QD nanostructure modelled 
by rectangular potential wells and a barrier. The barrier thickness was fixed as ds = 10 nm, 
the top barrier energy was changed according to Eq. (JHJ), the depths of both the potential 
wells were equal (Vl — Vr — —200 meV), and the thicknesses of both the potential wells 
were also the same di = d^ = d. We change potential well width d in the interval d e [6, 50] 
nm. The size of the region, in which the spin exchange appeared, has the width 2d + ds- 
Figure 7 depicts the calculated time duration of the first and second ix pulse as a function 
of width d of the potential wells. We see that the time duration of spin interchange quickly 
increases with increasing QD size. The n pulse time changes with the QD size as follows: 
T 71 ~ d 3 for d < 20 nm and T n ~ d 6 for d > 20 nm. These results show that even for the 
QD with the size such small as d ~ 50 nm the full spin swapping can be not realized during 
the coherence time. 

In the Heisenberg model of electron-electron interaction, the rate of spin interchange is 
determined by the exchange coupling constant, which is defined as the energy difference 
between the triplet and singlet states, i.e., J = Ex — E$- For the two identical QDs 
separated by the sufficiently thick and high potential barrier the triplet and singlet states 
are degenerate and J = 0. If the potential barrier becomes thinner and lower, the wave 
functions of the electrons localized in both the QDs start to overlap and the energies of 
triplet and singlet states become different, i.e., J ^ 0. The coupling constant J can also be 
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FIG. 7: Time duration T w of the first (n = I) and second (n = 2) it pulse as a function of thickness 
d of the potential well. Inset zooms in the short time part of the figure. 



n 

changed if the QD potential wells are different, i.e., the QDs are asymmetric |30j] . The larger 
the asymmetry of the potential wells the larger coupling constant J. Therefore, one could 
expect that the asymmetry of the QDs should lead to a shortening of the spin interchange 
time. However, it turns our 31[ that this is not always the case. 

We have investigated the influence of the asymmetry of potential well thicknesses on the 
duration time of spin interchange. For the fixed thickness of the left QD (di = 10 nm) 
we have performed the calculations for several values of thickness d^ of the right QD. The 
results are plotted in Fig. 8 for d^ = 5 nm and d^ = 15 nm. For sake of clarity, we present 
in Fig. 8 only the results for the electron in the right QD. The spin expectation values (Sf) 
for the electron in the left QD can be obtained by the reflection of S^(T) curves with respect 
to the axis = 0. 

The results of calculations are displayed in Fig. 8. If the right potential well is thinner 
than the left one, energy difference AEfi between the final and initial state tends to zero 
for long duration time T of spin interchange [dashed curves in Fig. 8(a)]. However, if the 
right QD is thicker than the left one, AEfi does not reach zero even for long time T [Fig. 
8(b)]. Moreover, it takes on negative values for some times T. The negative values of AEfi 
can be explained by the shift of electron density toward the wider QD in which the final 
state energy is lower. Analyzing the plots of Sf vs T we have observed that the accuracy 
of the interchange of spins reaches 99.9 %, but is never equal to 100 %. We conclude that 
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the asymmetry of the QD confining potential results in the increasing multiplicity of spin 
exchange, which leads to the longer duration time of spin interchange. 

(a) (b) 




T [ps] T [ps] 

FIG. 8: Expectation value = of z spin component detected in the right QD at the end 
of the process of switching on/off the exchange interaction (solid curves) and energy difference 
AEn between the final and initial states (dashed curves) as functions of duration time T of spin 
interchange for asymmetric QDs. The thickness of the left QD is fixed di = 10 nm, while the 
thickness of the right QD cIr = 5 nm in (a) and d,R = 15 nm in (b). 



V. LATERALLY COUPLED QUANTUM DOTS 

The electrostatic QDs with lateral interdot coupling are the subject of many recent studies 



0, [28|, [29|. In the laterally coupled QDs, the profile of the confinement potential can be 
modified by varying the voltages applied to the gates. This provides a fairly convenient way 
of tuning the potential barriers and wells. However, due to technological limitations, the 



size of the lateral QD is of the order of 100 nm [ijj], i.e., is larger than the size (~ 10 nm) of 
the vertical QD 27| measured in the growth direction. The nanostructure consisting of two 
laterally coupled QDs separated by the potential barrier has a typical linear size 200 — 300 
nm. Due to this relatively large size, the duration time of the spin rotation is expected to 
be long (cf. see Subsection IV. B). 

We still apply the theoretical model described in Section II with the frozen transverse 
motion of electrons. Now (y, z) are the transverse coordinates and we assume that the 
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electrons move in the x direction, e.g., the x axis can be directed along the effectively 
one-dimensional flow of electrons between the two lateral QDs jl4j ]. The larger size of the 
lateral QDs in comparison to that of vertical QDs leads to the weaker electron localization 



and the weaker electron-electron interaction. According to the theory [2l| the strength of 
the potential confining the electrons in the transverse directions determines the strength 
of the effective electron-electron repulsion [Eq. (pQ)]. [We note that now U e ff is a function 
of x.] The excitation energy hu± is a measure of the strength of the transverse parabolic 
confinement. In the calculations for the lateral QDs, we take on hu± = 5 meV. Then, the 
transverse confinement potential rather well approximates the confinement in the double 
lateral QD fabricated by the TU Delft laboratory [33|. The present model is also applicable 
to quantum wires, which are composed from different semiconductors, e.g., InP and InAs 



341 ]. that form a structure of potential wells and barriers. In this case, we deal with quasi- 



one-dimensional coupled QDs. In the present section, we modify the exchange interaction 
between the electrons localized in laterally coupled QDs by tuning the potential barrier as 
well as the potential wells. 



A. Symmetric quantum dots 

The profile of the confinement potential for the laterally coupled symmetric QDs is as- 
sumed in the form (Fig. 9) 

V(x) = kx 2 + V 2 exp [-(x - x ) 2 /d 2 ] (15) 

where k determines the strength of the parabolic confinement potential, parameter V 2 is the 
potential barrier height for V 2 > (potential well depth for V 2 < 0) , xq is the position of 
the center of the potential barrier (well), and d is the range of the Gaussian potential. 

The potential barrier separating the two QDs can be modified by changing V 2 [Eq. (fT5|) ]. 
In the calculations, parameter V 2 was changed in time in a smooth manner, i.e., according 
to Eq. (IHj), from 9 meV to and next from to 9 meV. Figure 9 shows the corresponding 
profiles of the confinement potential. We see that when changing the potential barrier we 
simultaneously change the potential wells. If the potential barrier is lowered, the electrons 
become localized closer to each other. The effective QD size (d^ R ) can be defined as the 
width of the corresponding potential well determined at the half of barrier height. At the 
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initial time moment, the effective size of the QDs df = df = 100 nm, the effective width 
of the barrier d R = 60 nm, the electron in the left (right) QD has spin +K/2 (-K/2). 




-200 -100 100 200 

z [nm] 



FIG. 9: Profiles of confinement potential E = V vs z = x for two laterally coupled QDs for V2 > 0. 
Parameters of the confinement potential [Eq. (|15p ] are k = 9.0 x 10~ 7 meV(nm) -2 , d = 185 nm, 
xq = 0, V2 = 9.0,4.3, and meV (solid, dashed, and dotted curves). 



In Fig. 10, we display the time evolution of spin indices S L,R and parameter V2. These 
results show that for the first 7r pulse with time duration Tf = 492 ps the spin interchange 
occurs with 100% accuracy. We see that during the larger part of time evolution the spins 
do not change. The spins start to rotate when the barrier height V2 decreases to ~ 5 meV. 
This suggests that the lowering of the potential barrier should lead to a shorter spin rotation 
time. If, however, the initial barrier height decreases, the effective QD size decreases too. 
For V2 — 5 meV the effective size of the potential wells and the barrier decrease to d e [ R = 76 
nm and df = 42 nm. This means that the shortening of the spin rotation time results from 
the size effect, discussed in Subsection IV.B. The fastest rate of spin interchange is observed 
for V2 < 1.5 meV. Then, the effective sizes are d e [ R = 40 nm and d% = 42 nm. 

If the QD size is of the order of several ten nanometers, the electrons are spatially sep- 



arated even in the absence of potential barrier [32] . For the QDs of this size the Coulomb 
repulsion is sufficiently strong so that the electrons are localized the QD boundaries and 
form Wigner molecules {32]]. Only the localization of electrons close to each other can lead 
to the rotation of their spins. In the model of nanostructure (Fig. 9), the effective width 
of the potential barrier is much smaller than that of the double potential well. Therefore, 
we can force the electrons to be localized in the same region of the nanostructure, if we 
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FIG. 10: Spin indices S p = S R for the right QD (solid curve) and S L for the left QD (dashed 
curve), and potential barrier height V% (dotted curve) as functions of time t for duration time 
Tf = 492 ps of the first vr pulse. 



allow the parameter V% in Eq. (|T5l) to take on negative values during the time evolution. As 
a result, the potential barrier converts into the potential well (Fig. 11) and the electrons 
become localized in the same QD for some time during the ir pulse duration time. 
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FIG. 11: Profile of confinement potential E = V [Eq. (I15p ] as a function of coordinate z = x for 
Vi > and Vi < 0. Shown are the plots for V2 = 9, 6, 3, 0, —3, —6, —9 meV (from top to bottom). 
At the initial time instant (solid curve), V2 = 9 meV, ds = 60 nm, di = dn = 100 nm. 

In the first computer run, parameter V2 was changed in a smooth manner from 9 meV 
to —3 meV. As expected, we obtained the considerably shorter duration time of the spin 
rotation: the first 7T pulse lasted merely for 46.1 ps, i.e., it was more than ten times shorter 
than in the case of V2 > 0. It appears a question: to what extent can we shorten the 
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duration time of spin rotation by lowering the potential well bottom? In order to answer 
this question, we have performed a series of simulations, in which the minimum value V 2 mm 
of parameter Vi was lowered from to —9 meV. In each simulation, the initial value of 
parameter V-i was the same: V%{t = 0) = V™ ax = 9 meV. The results are shown in Fig. 12. 
According to our expectation, the spin interchange duration time becomes shorter if V^ mm 
decreases. However, we also observe an undesirable effect, namely, the process ceases to 
be adiabatic, which results from the increasing amplitude of changes of V2 during time T. 
The non-adiabaticity is visible in Figs. 12(a-d) as oscillations of energy difference AEfi(T) 
and deformations of plots S^(T) for small duration time T. The process can be treated as 
adiabatic for T > 25 ps (T > 55 ps) if V^ in = (V^ in = -9 meV). We conclude that 
the lowering of V™ in below —3 meV is not advantageous. The shortest time of the full 
interchange of spins has been obtained for V™ m = —3 meV. Then, the orientation of spins 
was changed in the single swap after time T± = 46.1 ps with accuracy 100 %. 

We note that for the laterally coupled QDs the change of the interdot potential barrier 
into the deep potential well means that a section of quantum wire is created below the gate 
electrodes, which are responsible for the formation of electrostatic QDs 33J. However, in this 
process, the spin qubits can be destroyed if the electrons will tunnel from the two-dimensional 
electron gas [33j to the QD. Therefore, it seems that the mechanism of conversion of the 
interdot potential barrier into the potential well will not lead to the required rotation of 
spins in the laterally coupled QDs. However, this mechanism can be applicable to the spin 
qubits in the vertically coupled QDs [26| and the coupled QDs formed in quantum wires 
34 ] . 



B. Asymmetric quantum dots 

The bias voltage applied to the coupled QDs leads to the asymmetry of the confinement 
potential profile, i.e., the potential well depths become different. This suggests another 
mechanism of swapping the spins, which can be modelled by raising and lowering the bottom 
of one potential well. Raising the potential well bottom of one QD will stimulate the flow 
of the electron to the second QD. Therefore, both the electrons will be localized in the 
same QD. The subsequent lowering of the potential well bottom back to the initial level will 
cause the flow of one of the electrons to the first QD. If the duration time of changing the 
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FIG. 12: Expectation value = of z spin component in the right QD (solid curve) and energy 
difference AEfi (dashed curve) between the final and initial states as a function of time duration 
T of the process, (a) V 2 min = 0, (b) V 2 min = -3 meV, (c) V 2 min = -6 meV, (d) V 2 min = -9 meV. 



potential well will be suitably adjusted, the electrons will swap their spins. This mechanism 
is advantageous since the electrons meet in the same QD. According to the previous results, 
the stronger electron localization leads to the shorter time of spin swapping. 

In order to simulate this process, we model the confinement potential by the linear com- 
bination of two Gaussians, i.e., 

V(x) = V 1 exp [-(x - xrf/dj] + V 2 exp [-{x - x 2 ) 2 /dj] , (16) 

where V\ and V 2 are the potential well depths, x 1 and x 2 are the positions of the centers of 
the QDs, and d\ and d 2 determine the sizes of the QDs (Fig. 13). 
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FIG. 13: Profile of confinement potential E = V [Eq. (I16p ] as a function of coordinate z = x for 
several time instants. Parameter V\ = —10 meV is fixed. At the initial moment (solid curve), 
V2(t = 0) = —10 meV. Shown are also the potential profiles for V2 = —7, —5, —3, and meV. 

We change parameter V2 according to the smooth function of time [Eq. (JHJ)]. At the 
initial time, the electron in the right (left) QD has spin +h/2 (—h/2). We have performed 
several simulations for different tc pulse duration time T w with the fixed parameters of 
confinement potential [Eq. ffTB"]) ] V\ = —10 meV, d\ = d,2 = 64 nm, and X2 = —X\ = 58 
nm. For t = (solid curve in Fig. 13) the bottom of both the potential wells was at 
the level —10.4 meV, effective widths d\ = d,2 = 48 nm, ds = 51 nm, and the relative 
height of the barrier was 1.7 meV. In each computer run, the initial (minimum) value of 
parameter V2 was V^(t = 0) = V 2 mm = — 10 meV. The potential profiles shown in Fig. 13 
suggest that it is not necessary to remove the right potential well in order to localize both 
the electrons in the left QD. Again, we are looking for the answer to question: for which 
range of changes of the right potential well bottom the duration time of interchanging the 
spins is the shortest? The results of the simulations are displayed in Fig. 14. As in the 
previous figures, for sake of clarity, we present only the spin expectation values for the 
right QD, since Sj? can be obtained by the reflection with respect to zero ordinate axis. If 
parameter V2 is changed from V 2 mm = — 10 meV to V 2 max = 0, the changes of the amplitude 
of the right potential well bottom are so large that even for T ~ 100 ps the spin interchange 
process is non-adiabatic. In this case, we observe the spin exchange with accuracy 99.8 % 
for T[ = 90.4 ps. The adiabatic process of full interchange of spins is observed for T of 
the order of several hundred picoseconds. Therefore, the process, during which one of the 
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potential wells is entirely removed, is not advantageous. Therefore, we reduce the amplitude 
of changes of parameter V 2 by decreasing its maximum value V™ ax . The useful results are 
obtained when raising the potential well bottom to V 2 max = —3 meV [Fig. 14(a)]. However, 
we have to take into account that V™ ax can not be too small since then the potential well 
bottom is not raised high enough so that the probability of electron tunneling from the right 
to left QD is too small, which leads to the long spin swap duration time. For example, for 

ymax = _ ? me y [p ig> 14 ^] yw = 71 g ^ shortegt ^ pulge t j me (rpn = 13 g pg ) hag 

been obtained for V™ ax = -5.5 meV [cf. Fig. 14(b)]. 
[Fig. 14(b)]. 

We conclude that the time changes of the confinement potential, which lead to asymmetry 
of QDs, result in the considerable shortening of the spin swap time if the amplitude of 
potential well changes is of the order of ~ 5 meV. 

VI. CONCLUSIONS AND SUMMARY 

In the present paper, we have proposed a method for a theoretical quantitative description 
of spin swapping process in a double QD. Using this method we have estimated the spin 
swap duration time induced by the exchange interaction between electrons in vertically and 
laterally coupled QDs. The results obtained are also valid for quantum wires, which contain 
coupled QDs. The switching of the exchange interaction, considered in the present paper, 
was triggered by the corresponding time changes of the confinement potential profile. The 
time dependence of the confinement potential is of crucial importance in the process of 
spin interchange. We have shown that the confinement potential should change in time in 
a smooth manner and sufficiently slowly to ensure the adiabaticity of the spin swapping 
process. 

We have also shown that the full spin interchange (with 100% accuracy) can hardly be 
achieved in the single spin swap process. In the coupled QD system, the full spin inter- 
change occurs as a result of many-fold incomplete spin swaps. The most promising results, 
i.e., the shortest time of full spin swapping in the adiabatic process, have been obtained 
when smoothly changing in time the shape of the potential in the interdot region from the 
potential barrier to the potential well. However, it seems to be difficult to realize this pro- 
cess experimentally in the same nanodevice by changing the voltages applied to the gate 
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FIG. 14: Expectation value = (solid curve) of z spin component of the electron in the right 
QD and energy difference AEjj (dotted curve) as functions of duration time T of potential changes 
for asymmetric QDs. (a) V^ ax = -3 meV, (b) V^ ax = -5 meV, and (c) V^ ax = -7 meV. 



electrodes, which define the laterally coupled QDs. 

The results obtained for asymmetric QDs are rather disappointing: it turns out that the 
asymmetry of the confinement potential leads to the increase of the duration time of spin 
interchange, since in this case the full spin swapping is obtained after several incomplete 
spin swaps. Both in the vertically and laterally coupled QDs the spin swap duration time 
decreases if the range of the confinement potential becomes shorter. This leads to the 
conclusion that the size of of coupled QD system should be smaller than ~ 100 nm in order 
to obtain the sufficiently short spin swap time. 
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In summary, the exchange-interaction induced processes of spin rotation can serve as 
a mechanism to manipulate with electron spin qubits in coupled QDs. The spin swap 
duration time, which is several orders of magnitude shorter than the spin coherence time, 
can be obtained for sufficiently small coupled QDs, in which the gate voltages cause the 
smooth time changes of the confinement potential. A fabrication of the nanodevices, which 
satisfy these requirements, is a challenging task for the nanotechnology. 
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